General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



; cf-m 7 77 

(NASA-CR-171777) DYNAMICS ANALYSIS OP N84-22598 

ELECTRO DYNAMIC SATELLITE TETHERS. EQUATIONS 

OF MOTION AND NUMERICAL SOLUTION ALGORITHMS 

FOR THE TETHER Final Report (Texas Univ. ) Unclas 

46 p HC A33/MF A3 1 CSCL 22B G3/15 19090 


Dynamics Analysis of Electrodynamic 
Satellite Tethers 


NASA - JSC Contract 
NAS 9 - 16744 

Final Report 
January 31, 1984 


"Equations of Motion and Numerical 
Solution Algorithms for the Tether" 


by 

Dr. Paul E. Nacozy 
Professor 

Center for Space Research 
Dept, of Aerospace Engineering 
and Engineering Mechanics 
University of Texas at Austin 
Austin, Texas 78712 

512 - 471-4248 



1 


Introduction 


In this report the equations of motion are developed for a perfectly 
flexible, inelastic tether with a satellite at its extremity. The tether 
i3 attached to a space vehicle in orbit. The techer is allowed to possess 
electrical conductivity. Also in this report a numerical solution algorithm 
to provide the motion of the tether and satellite system is presented. 

This report has several purposes. 1) To provide a solution algorithm 
for the motion; 2) To allow an analysis into the physical and dynamical 
properties of “the motion for various electrical currents and other 
parameters; 3) To provide a check on the developments and algorithms of 
other researchers; 4) To allow a determination of the approximations, if 
any, introduced in the developments of others; and, 5) To provide a set 
of exact, analytical differential equations that describe the motion. 

The purpose of having exact, analytical differential equations has 
se>eral benefits. The resulting differential equations can be solved by 
various existing standard numerical integration computer programs. The 
resulting differential equations allow the introduction of approximations 
that can lead to analytical, approximate general solutions. The 
differential equations allow more dynamical insight of the motion. 
Modifications, extensions, generalizations of the dynamical system are much 
easier to accomplish for the system of exact differential equations. And 
many numerical integration errors are more conveniently monitored and 
controlled. And finally, overall comprehension of the equations and 
solution algorithms is somewhat more straightforward. 
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2 . Equations of Motion for the Tether- Sa tellite System 

Consider an element or a segment of the tether with length A S and 
mass per unit length along the tether, p. Consider a tether having a 
negligible diameter relative to the length of the tether. 

Throughout this report, the tether will be considered as inextensible 
or inelastic. 

Figure 1 shows a portion of the tether, separated into several, 
adjacent, elements. 



The force of one element on an adjacent element will be denoted by the 
vector F. Since the forces between two adjacent elements (at the junction 
point) are equal e.nd opposite, F on one element (1) due to an adjacent 
element (2), gives -Tf for the force on (2) due to (1). Let s denote the 
distance along the tether, where s - 0 at the point of attachment at the 
space vehicle and s - L at the satellite, where L denotes the length of the 
tether. Then 1" will be a function of s p R « R (s). The function R(s) will 
increase or decrease continuously as a function of s from one element to 
the next, as illustrated in Figure 2. 
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Consider the forces acting on an element as shown in Figure 2. 



0 


Figure 2 

The forces and -Rg are the forces due to each adjacent element. 

ft ft 

The forces k (s^, k (s 2 ), are other external forces such as 

gravitational forces, electrodynamical forces, etc., acting at various 
locations s i» s 2 » s 3 » .... on the element where s is the distance along the 
tether. 

ft ft ft 

Define vectors r., r D , and r , from a reference origin 0 to 
fl o cm 

each end of the element and to the center of mass of the element, cm, as is 
depicted in Figure 2. Note that the point cm is not necessarily situated 
on the element. 
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The first equations of motion for the element A S is the sum of the 
forces acting on A S equal to the time rate of change of the linear angular 
momentum of the center of mass of the element, cm, which, in this caso, is 
equivalent to the mass of the element times the acceleration of cm, or 


A S 


..* 


cm 



n 

+ I 
i-1 


0 ) 


where dots denote differentiation with respect to time, and p is the mass 
per unit length of the tether. 


The restriction placed on the origin 0 is that it be an inertial 
point. 


The second equation of motion is the sum of the moments about cm 
equal to the time rate of change of the angular momentum about cm, or 



N 

Z 

i-1 


flcm 


( 2 ) 


The forces k (s) in equation (1) are not discrete forces but are 
distributed uniformly and continuously over the element, hence the 

_# — 

summation given in equation (1) over k (s) can be replaced by k(s), and the 
total force is determined by. 


k(s)ds 
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where the limits indicate that the integration is taken from r^to r g and 

cm 

similarly the moments due to k (s) about the point cm, M A in equation (2), 
are: 


!a ( r s ' r on,> x k(s)as 
_# 

where r 3 is the vector to the point s on the tether. 

Equation (1) becomes 

..* 

1 S • 0 ■ F on. ■ \ - >>B * <3) 

and equation (2) becomes 

& * 3 - c c 3 xi b * f» \ as * is as <“) 

where A is the angular momentum of the element at s, with length A S, 
about the center of mass of the element. The quantities and M g are the 
torsional moment and the bending moment, respectively, per unit length. 

The vectors R A , and -R g can each be decomposed into two component 
vectors at right angles to each other, one tangent to the element and one 
perpendicular or normal to the element, as shown in Figure 3. 
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The vectors T A and -Tg, components of R A , and -R g along the element are 
referred to as tensile stress. The vectors Q. and -Q q are the shear 
stress, since they are normal to the element. 

Dividing equation (3) by A S and using the component vectors of 
and Rg gives 



Letting A T - T a - T 0 and A Q - - Qg, and taking the limit gives 


11m AT AQ JT 3Q 
AS -*> 0 "AS " AS “ 3s ” 3s 

Also 



where K(s) is a force per unit length at s. This result follows since as 
A S-*0 , this implies that point A (or point B) approaches, in the limit, 
point B (or point A). 

Equation (5) becomes 


• • 


« 


P • 


r 


cm 


3T 9Q 

as as 


♦ K(s) 


( 6 ) 


6 


ORIGINAL PAGE f9 

OF POOR QUALITY 


Equation (6) is the translational equation of motion (for the aenter of 
mass) of the inf lnltessimal element ds, Including torsional and bending 
moments, for a non-flexible, inelastic tether. 

In equation (4), dividing through by A s and U3ing the component 
vectors for W A and Rgi 


e (4I W (r A - r *„ ) , (r I - r =n, ) , r A - r cm 

3t~S 55 xT a * 55 x0 a 55— xT b 


_# _# 

" B " x5 B * 55[a <? ( S ) - F cm ) x * (s)<ls 


‘ 55 fS V 3 * 5$ V’ 


(7) 


Taking the limit of equation (7), as A s-»0 gives: 


lim 


A S-*0 dt 


A H 


cm 


A S 


d cm 
dt 3S 


-* _ -* 

— ilE_ r , A —- c J}. xT . dr_ t 
AS-K) AS x A ds A u 


* « 

lim B cm *t _ dr ™ * 

AS'-^O AS B ds B u 


7 


QKHSiMAL PAGE 
OF POOR QUALITY 


and letting Q A and Q 0 go to Q as A s 0; 


11m 

* 

ft ft 

r A " r cm 

AS-*0 

AS 


* « 

r A ' r cm 


XQ A + AS 


X(-Q B ) 


11m 


AS-*Q 


ft ft « ft 

"jLlIom I q« J_ r B ^ 

AS W A AS XW B 


11m 

AS*»0 


—ft _» _« —* 

r B ' p ora r om • r A * 

Is xQ 


11m 

AS-*«0 


ft ft 

r D - r 


B —1 A - d/ X Q 

AS w ds w 


where the direction r 


A 


—» 
to r B 


is 


chosen as the direction of 



11m 1 f B , » -* . 

SO AS)A ( r(s) " r cm ) xK(s)da 


AS- 


since r(s) - r - 0 and ds/AS - 1 in the limit. Also, 


11m _1_ 

AS^'O AS 


H f 


A M T 113 • «T 


and 


11m J_ 
A$-*0 AS 


\l% 


ds - M B 


So equation (9) becomes 


Also, 


0 
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d 

3H-_ 

cm 

dt 

3s 


Equation (8) 
element ds 
moments, for 


• - *5 * m t * m b < 8 > 

is the rotational equation of motion for the infinitessimal 
about its center of mass, including torsional and bending 
a non-flexible, inelastic tether. 
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jU jl Perfectly Flexible Tether 


The torsional moment M T , Is the resistance of the element to twisting 
about Its axial direction. In this report, Inclusion of the torsional 
moment Is not necessary for the tether. Hence ■ 0 in equation (8). 

The bending moment Is the resistance of the element to bending about 
an axis perpendicular or normal to the physical element. The bending 
moment Is determined by the properties of the tether. For this report, it 
will be assumed that the tether is perfectly flexible, or has no stiffness, 
and that the bending moment is negligible. This is probably a good 
assumption for many types of long, thin tethers. However, if the diameter 
of the tether were appreciable, this assumption might need evaluation, 
depending on the physical properties of th<s tether. Hence, in this report, 

K p - o 

in equation (8). 

For a perfectly flexible tether, with no bending moment, the vector Q 
in equations (6) and (8) will be set to zero. (See, for example, Shames, 
Engineering Mechanics , Statics and Dynamics , Prentice Hall, 2nd edition, p. 
148, for the proof that O’ - 0 for a perfectly flexible cable in equilibrium 
(zero acceleration). Extension of the proof to the case at hand, the non- 
equilibrium, non-zero acceleration case, has not been given. However, it 
will be assumed here that for small to moderate accelerations of the 
tether, setting Q - 0 will provide sufficient accuracy.) 
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P r 


cm 


3T 

TG» 


K(s) 


( 9 ) 


j 3H 
d cm 

it TT 


( 10 ) 


In equation (9). the vector T lies along the tangent to the tether at 
the point s, as discussed above. Equation (10) can be integrated 
immediately giving 


» constant » c(s) (11) 

providing the result that the angular momentum of each (infinitesimal) 
element of the tether remains constant in time. Note that the constant in 
equation (11) is only "constant" in time, but in general is a function of 

s; the value of the "constant" C(s) depends on the initial conditions. Of 

course the reason that we have this result is that we set Q - M_ * H r - 0, 

o T 

which are all the moments about the center of mass of ds. 

Equation (9) is the translational equation of motion for each 

infinitesimal element, ds, for a perfectly flexible, Inelastic tether. In 

m g ^ 

equation (9), replace r by r(s.t) ,T by T(s,t)and T?(s) by K(s,t). And 

era 

to denote that the total derivatives with respect to time are to be taken 
with s held constant, the total derivatives are replaced by partial 
derivatives. Equation (9) becomes 


- - 3 T ( , 3 . ! - U ♦ K(s,t) 


Dt 


Ds 


( 12 ) 


where K(a,t) represents the external forces per unit distance along the tether 
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!• The Inelasticity constraint 


The assumption that the tether Is inelastic (no stretching) can De 
written as 


3r 


3s 


• 1 or 


£•£ • - 


(13) 


Moreover, the vector Ar AS, in the limit as AS goes to zero, along with 
constraint ( 13 ), coincides with the unit vector along the tangent of the 

A 

tether at the point s. Denoting this unit vector by t 




Equation (12) then becomes 


(14) 


p r(s,t) # - - 
3t Z 


■3s 


T(s, t) 


3~(3,t) 


as 


♦ K(s, t) 


(15) 


The inelasticity contraint provides a differential equation for the 

3r* 

determination of yj-and T. Differentiating equation (13) twice with 
respect to time gives 


3r 


_» 
3 3r 


3s * 3T IT 


( 16 ) 


and 
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3 

3r 

3 

3r 

3r* 3 2 

a» * 

... * l 

3t 

3s 

* 3t 

3s 

33 * 3t 2 

3sJ 


( 17 ) 


_# 


Also since r (s,t) and its derivatives are assumed to be continuous, then 


3 

_* 

3r 

3 

_* 

3r 

3t 

3s 

" 3s 

3t 


( 18 ) 


and Equation (17) can be written as 


3s 


_# 

In 

3t 


_3_ 
' 3s 


_* 

3r _3_ 
3s * 3s 


2 -* 
3 r 


3t 


(19) 


Substituting equation (15) into equation (19) gives; 



- - 

* 

— 

_# 

_* 

3 

3? 

3 

3r 

1 3r 3 r _ JL 

p 3s * 3s L 3s 


3s 

3t 

‘ 3s 

3t 

— 

* ^ 

OJ 


+ K] - 0 


or 








— — 



3 

3r 

3 

_* 

3r 

_* 
1 3r 

jL 

T £- 

_# 

+ liL. 

3K 

3s 

3s 

’ 3s 

3t 

p 3s 

3s 1 - 

T 3s 

p 3s 

3s 


( 20 ) 
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5_. The Equations of Motion for the Satellite and the Space Vehicle 

In this section, differential equations will be developed for the 
satellite and the Space Vehicle in order to determine the motion of the 
entire satellite - tether - vehicle system as well as to provide 
information concerning the boundary conditions for the tension at each end 
of the tether. 

The forces acting on the satellite are the tension due to the tether 

T(l,t); and other external forces, 1^(L,t), where L denotes the value for 

s at the end of the tether, ard K^L.t), denote forces acting on a point 
mass as opposed to K(s,t) which is a force per length. 

Hence , 


?— * 

3 r(l,t) 

3t 2 


m 


-1 

SAT 


T(L.t) 


3r(L,t) 

3s 


lv( L, t) 


SAT 


(21 ) 


or 


p * 

3 r(L,t) 


3t 


T(L,t) ♦ K^L.t) 

m SAT 


where, as discussed in section 4, 


( 22 ) 


T(L,t) 


A 

T(L,t) & ■ 


_* 


(23) 


and where the negative sign on the first term on the right in equation (22) 
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indicates that the force of tension on the satellite is in the negative s 
direction, toward the tether. 

For the other end of the tether , attached to the space vehicle, make 
the following assumptions. Assume the tether is attached to the space 
vehicle at some point 0*. The point 0* is in general not an inertial 
point. Let 0* be the origin of our coordinate system with rectangular 
Cartesian coordinate axes defined as follows. The z-axis is along the line 
from the center of the earth to the point O', positive upwards. The x-axis 
is in a direction such that the x - z plane contains the velocity vector of 
the point O', with positive direction in the direction of the velocity 
vector. And y-axis is defined by the right-hand rule. 

This definition of the x-axis is convenient, since it is in the 

direction of (but not necessarily coincident with - allowing for non-zero 

eccentricities of the space vehicle’s orbit) the velocity vector of the 

point of attachment at the space vehicle. For zero-eccentricities and 
assuming O' at the center of the space vehicle, the velocity vector will, 
of course, lie along the x-axis. And the y-axis (determined by right-hand 
rule) will be normal to the orbital plane. For the actual situation, 
assuming the space vehicle is rotating at a reasonably slow rate, and for a 
non-zero eccentricity, and including perturbations, the x-axis should be 
close to the plane containing the position and velocity of the space 
vehicle. Let the unit vectors along the x,y,z axis be !,J, and k, 

respectively. 

The coordinate axes will be rotating with angular velocity a? . The 

* 

point 0* will be located at s - 0 or at r(0,t) . The velocity and 
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acceleration of the point of attachment will be denoted by 


3r ( 0 , t ) 
3t 


and 


a 2 r(0,t)* 

2 1 
3t 2 


_ * 

respectively. Let r (t) be the vector from the center of the earth to 

S • V • 

a reference point within the space vehicle. The acceleration of the 
vehicle is 


52r 3.V. (t) 


2t 


m 


— T(O.t) 


s. v 


acting at 0* 


+ 


1 


m 


s. v. 


K 


s. V. 


(t) 


(24) 


where T(0 f t) is the force on the vehicle due to the tether at s - 0 or at 
O' and K_ .. (t) are all other external forces such as gravity, drag, 

o • V • 1 

_ ft 

thrusting, etc. acting on the space vehicle, and located at r (t) . 

S • V • 

The acceleration at the point of attachment, O', is 


3 2 r(0,t)* 


» 2 r 3 . v . (t)1 


s2 r rel (t) 


3t 


3t 


3t 


(25.1) 


where 


t. r rel (t) 
3t 2 


“s.v. x(r<0 ' t >* - r 3 .v. (t),) 
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i.v. 


x(r(0,t) 


* — 


- r 3.V (t>#)) 


( 25 . 2 ) 


where ft and Q are the angular velocity and angular acceltration of the 
space vehicle, respectively. 


In some cases, for approximate solutions, some assumptions can be made to 
decrease the time of computation. For example, assume that m is very 

3 • V « 

large compared to the tether. Let this be assumption A. And assume that 


s. v. 


T(0, t) « 


1 ^ 


3. V. 


3.V. 


(t). 


Let this be assumption B. 


And if one assumes two-body motion for the space vehicle’s center of mass 
(assumption C), then equation (24) becomes 


3t 2 


-Gm earth r s.v. (t) 


r s.v. (t) 


( 26 ) 


The assumptions that the mass of the tether and that T(0,t) on the space 
vehicle can be neglected for the motion of the space vehicle (assumptions 
A+B) allow us to be able to determine a priori the acceleration of the 
space vehicle,, and, as will be seen later, allows us to determine T(0,t) 
readily, for the tether. 


The assumption that IC( 0 , t ) contains only two-body forces (assumption 
C) allows us to obtain the solution for the motion of the space vehicle 
analyticaly, for two-body motion. Hence, 

r(0,t)\ |f (0,t)*' -2— (0.t>\ 

at 2 


as well as u(t), oj(t), are known a priori, with these assumptions . 
Another assumption can be made by placing the point of attachment O', 


IS 


—ft 

at the center of the space vehicle, at r (assumption D). If the errors 

s • v • 

made by the above assumptions are larger than desirable, then the 

g 2 r ft 

determination of (Q,t) will not be known a priori and the motion of 
the satellite - tether - space vehicle system must be solved for 
simultaneously (see for example report #1 of this contract). The 
developments, equations and algorithms given in this report do not contain 
any of the assumptions just discussed, but are for the general situation, 
where the motion of the space vehicle and the tether are solved for 
simultaneously and all external forces on the tether and space vehicle are 
included. 

_ * 

Let r(s,t) be from 0 to the point s on the tether, then 


r(s, t) 


r*(s, t) 


r 


(s-0, t) 


(27) 


where r (s - 0,t) is to the point of attachment O', since s»0 at that 
point. Hence r(s,t) is the vector from the O', at the vehicle to the point 
s on the tether. And, for a non-inertial, rotating reference frame: 


9r(s,t) 


ft 


3t 


9r(s,t) 

2t 


+ a x r(s,t) 


(28) 
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3 2 r(3,t)* _ a 2 F(0,t)* + 3 2 F(s,t) 


3t' 


at' 


at' 


♦ to x( uT x F(s,t)) ♦ to x r( s, t) 


♦ 2 to x 


3r(s,t) 

3t 


(29) 


The first term on the right can be given by equations (25.1) and 
(25.2), the fourth term would be included for non-zero eccentricity orbits 
of the space vehicle, the third term is the centripetal acceleration, and 

a 2 pV a t) 

the last is the coriolis acceleration. The second term, • •• • • - * - is the 

3s 2 

acceleration of the element at s with respect to the point of attachment , 
£' . Notice that the last term contains a velocity with respect to the 

point of attachment, O', that is, the relative velocity, defined 

in equation (28). 

If an additional assumption of zero eccentricity is imposed 
(assumption E) along with assumptions A, B, 4 C, and also with the point of 
attachment placed at the c.m. of the space vehicle ( assumption D) then : 

F(0,t)* - constant . f? (30) 
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and 
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3r(0,t)* Inearth 


3t 


r(0,t) 


v circular * 


and 


— circular 

« * — 

r (0,t) 


co * 0 


and eqt. 25.1 becomes simply 


(3D 


(32) 


(33) 


(34) 


3t 2 at 2 

In general (with no assumptions) and with the reference frame defined 
above, equation (22) becomes, using equation (29) 


r(L,t) 

3t 2 


-1 

m SAT 


T(L,t) 



(L,t) 


2 — * 

3 r(O.t) 


at' 


- w x co x r(L,t) 


- w x r(L,t) 


2d) x 


3r(L,t) 

at 


(35) 
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with assumptions A,B,C,D,E, (neglecting tension on the space vehicle, 
placing O' at the center of mass of the space vehicle, and assuming 
circular motion reduces) equation (35) to 


2 

r( L,t) - ~ - T(L,t) + -J—K p (i,t) 
3t 2 m SAT m SAT 


Gm earth r(0,t) 

+ — -- - « x (o» x r(L, t) ) 


O(o, t )*] 3 


2 us x 


3r(L,t) 

3t 


(36) 


where r(L,t), , -—s- r (L,t) are with respect to the rotating 

at 2 

coordinate system with origin at the center of mass of the space vehicle. 

* 

Note that in equation (36), the term with Z does not appear since with the 
assumptions made, ^ - 0. 
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6. Equations of Motion In the Rotating Reference Frame 

In the rotating coordinate system, equations (12) (or (15)) and 
equation (20) must be modified. 

In equation (12), the acceleration is replaced using equation (29), 
The result Is, using T - T || ana * T & , 

6 3 


5 2 r(s,t) 
2t 2 


■ I 3T 3i* 1 T 3r 1 tt, , 

f'St ■ 35 P ^2 * f 


2 . » 

3 r(0,t) — 

to x (to xr(s,t)) 

at^ 


(37) 


- to xr(s,t) - 2 Z x |£(s,t) 

For equation (20), rederive equation (19) by taking partial 
derivatives of equation ( 13 ) with respect to the rotating reference frame 
centered on O'. The result is 


3r 

3 

3r 

. 3r 

a 


at 

* 3s 

at 
» » 

3s 

’ 3s 

1st 2 ] 


as 


Substituting equation (37) into this result gives 


3 _ar 5 3r 3r r 1 3 * . 1 rr 3 2 r J 

as at ‘5s at as as L p "5s T p K 3t 2 


(37,1) 


- u j x (w xr) -toxr-2tox|£]-0 


23 


ORIGINAL PAGE IS 
OF POOR QUALITY 


and taking the partial derivative of the terms within the brackets 
provides: 


3 , 3r v 3 

13 ( It } 13 


<£ > 


I in 

p 3s 



M 

3s 


3r 3 S 1 2 ? 

s5 7»^2 <0 - t > 


3r 

3s 


u x(b) X 


3r, 

13 } 



CD X 


3r _ 3r 
3s 3s 


2 


u x 


3 

3s 3t 


0 


The term . u x upon rearrangement, vanishes. Using T ■ T and 
* • 

3 — 

since r (0,t) - 0, gives 


3 

13 


3r 

3t 


3 

1? 



1 ifl 3r 
P 3s 2 38 


3r 2 3r 3^ 3T 
3s " J Is ‘ . 2 13 


1 T & 3r 1 3r 3K 

p as 3 ’ 3s p * 3s ‘ 3s 


Either equation 


3r — _ 3r . 3r 

" 13 * “ x * Is } "13 


(20) or (38) can be used. 


. 2 a) x 
However 


3 3r 
Is 3t 


0 


, if equation 


(38) 

( 20 ) 


dp 

is used, must be obtained using equation (28). Since relative 
velocities and accelerations will be used here, equation (38) is more 
convenient. 


For the satellite, the equation corresponding to equation (37) is more 

3 2 f 

complex since the quantity ■ s —r is difficult for the satellite and the 

3t 

d dp 

adjoining element of the tether. The terra A becomes 

38 3t 2 
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3 3^F(L,t) 3 "r(L.t) a^s.t) . n . 

K 9t 2 -* L 9t 2 9t 2 


Substituting this definition into equation (19). along with equations (35) 
and (37), gives 


3 ,3r\ 3 . 3r ^ , 1 -1 r n ¥S 

sT'ii 'it' ts • sc ( ts i-t ‘"• t ’ 

Ski 


rrpi^ (L.t) - « X G xr(2,t)) * 2 m x - « xF(L,t) 


mSTT 


♦ j T(a,t) - i lf(«,t) ♦ « x(w xr(s,t)) 


+ a x r(s,t) ♦ 2 « x ® 


25 


ORIGINAL PAGE IS 

OF POOR QUALITY 

rearranging terms gives 

J_ 3 r<£E\ . 3r lia , 1 v r 1 . 

7s ( f ) * -5s ( F ) * is • rrW c * s^: T(L » l) 

♦l^-Ks.t) ♦ r L_K*(L f t) -l?(s # t) 
e rs ™SAT p 

- « x (w x[r(L,t) - r(a,t)]) « 2 « x 


^ • 

C^ t) * $r(a,t)3t] - w x[r x(L,t) - r(s,t)] - 0 


Taking the limit, where feasible, gives 


3 ,3r 


11m 


r 3r.\ 3 / 3r x 5r * t 1 x r 1 ar#, 

W 'it 1 • -5? < jt > • Ul ( CT )C T(u 


t) 


SAT 


*7^ 1< ‘- t) * ITT 


SAT 


_ in _ 3r .-..a 3r(L, t) «r - ..3r(L,t) 

3a * “ x ' u * 5 T ' ~ • xrr — *7 rr w x — rr — 


3s 3s * ™ x 7s 3t 3s 


3s 
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The last term vanishes, and the final result is, 


3 ,3r N 3 ,3r x A 3r lim , 1 


Q ( 01 \ O / < 7 » \ . Ol / I \ r 

1$ { W •'§I ( 3t )+ 3i* TK ( t=a } • C * vu ’ w T TIE 


_1 T (L.t) ♦ T(a,t) 


♦ ML.t) - 1 K(s, t) ] - 3r( j; ,t) (L.t) . O) x(u x|£ (L.t) 

p dS d3 


3r* / 1 \ ^ « 3 3i* ( L 1 1 ) 

te <L ' r) • 3 “ *5? —jit 0 


or, 


3 / 3r ^ 3 5r A 3r lim^ 1 . r ~1 — / , *. \ 

"5? ( 3t } *?s3t to ‘ a*i?I=s' [ 3i — r(L,t) 


SAT 


+ • Ir + 7 T + 7~ — K*CL,t) - 1 K(s, t) 

p 3s 3s P 3s 2 at P 


3r( L, t) - „3r(L,t) 3r(L,t) „ - 3 3r(L,t) rt 

~T~ ‘ “ X ~3i 3s ' 2 w X 31 3t * 0 


(39) 


For the space vehicle, the corresponding quantity is: 


3 3 2 r(0,t) 
3s 3t 


lim 
s-> 0 


1_ 

s 


3 2 r(s,t) 


3t 


3 2 r(0,t) 


Dt e 
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Substituting this expression into equation 37.1 using equation (37) gives 


3 3? 3 3r 3r 

"5s * at * "Ss 3 t "5s 


* 

k 


* 

J K(s,t^‘- u x(u xr(s,t) - U) xr(s,t) 
« 


-z - &< ■ %» - . o 

at 2 at 2 


j2 r / Q 

ITie quantity is the relative acceleration of 0* at 0* and 

3t 


hence is zero. 


Use the following! 


AAIU . — . 

3 ^Q J [ w X ux F(s,t) 3 - U) X U x a - 0 ,t) 


since r(0,t)« 0, and similarly 


lira . , — 

Trt C “ x r(s * t >> «■» ^ ' 


lira . _ , . 

~ 2 . r 2 iii y ( s 1 1 ) i _ — a 

s->os L x at ■* 2u)X 'aF 


iliiffg fc ) 
3*0 S rlS,U 


■ 2 'x j- i r .( _o, t ) . 2 u .. a ar(Q.t) 
w at 3s 2 “ x si" at 
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Equation (40 becomes) 




||. uxux iE^,.|E. 2 ux ^£^t) 


-in -vMMl _ 3r 3r*(0,t) . 

33 33 * ’ ' 0 


as 


where the last term vanishes. 


(41) 
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Numerical Solution Algorithm 


A. Relations and Definitions 

Define the following quantities: 

r(s,tp • Position of element at s at time t i relative to space vehicle 
(O’) 


aris.tj) 

n ^ ■ Velocity of element at s at time t^ relative to space vehicle 

(O’) 


3r(s f t^ 

^ - Unit vector having a direction of the tangent at s at t^ 

a^s.tj) 3 3 r(s f t 1 ) - 

, r- First and second derivatives of -r- 

S3 2 33 3 33 

3 3r(s,t.) 3r(s,t.) - 

^ (*gg- ^ ) - rate of change of the unit vector or with 

respect to time or rate of change of the velocity vector with 
respect to s - all in the rotating coordinate system. 

K(s,t.) - all external forces on s exceot for tension, per unit length. 

_ # _ # 

K(0,tj) and k(l,t.) = external forces on the space vehicle and 

satellite, respectively 


A(s,t i ) 


1 3r(s,t.) 
7 3s 


( 46 ) 
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AA(s, 


BCs.tj 


CCs.t^ 


D(l,t. 


E(s,tj 


G(a,t i 


H(s,t i 
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1 d 2 r'(a,t l ) 

’i ] ' T 3 


3s‘ 


( 47 ) 


1 _ _ _ 3 r(s,t ) 

)- - Ms.tj) - u) x (w x r(s,t.)) - 2 u) x - u x r(s,t i )( 48 ) 


3 r(s,t, ) 3 r(s,t.) 

)- 4 ( iir-i- ) . 4 ( =r— !-) 


3 s 3 t 


3 s ' 3 t 


( 49 ) 


3 Srts.tj) 3 r(s,t i ) 


3 t 3 s * 3 t 3 s 


(50) 


^ * i -1 *- k(l,t.) - 2 u x - a) x(o) x r) - w x r(s,t.) ( 51 ) 


“SAT 


3 t 


)- 


T 3 k(s,t.) 3 r(s,t.) 


P 3 s 3 s 


3 r(s,t ) _ _ - 

— §i — • “ * <» x H> 


( 52 ) 


ar -i S - t i ) 2 - x 3 3f 
3 s * 2 “ X 3 s 3 t 


( 53 ) 


3r(s, t . ) 3r(s,t ) 

I = J i_ L_ 

p 3s 3s 


( 54 ) 


1 3 r(s,t i ) 3 3 r(s,t i ) 
“ ~ “ "as" 3 " ' 


( 55 ) 
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EE(L,t i )- - 


3r(L,tj) 

3s 


u x (u x 


3r(L,t 1 )) 

3s 


FF(s, t. )- 


3r(L,t i ) 

3s 


-2 3r 
p * 3s * 


. 2 u x 



3 3r(L,t) 
"5s 3t 


( 55 . 1 ) 


( 55 . 2 ) 




At time step t^ f assume the following are known from the previous 
integration step, or given as initial conditions and let L be the length 
of the tether: 

_ SrCs-.t.) 

r(s j ,t i ) ' §t 

for 0 S Sj S L;j * 1 ,2, . . . ,n 

The quantity Sj indicates the center of the j th element, as shown in 
figure 4. The integer j indicates the segment, beginning with s, adjacent 

to the space vehicle, with n elements. The vector F(Sj,t i ) is the 

fVom - 

position vector to s. the point of attachment, O'. The vector -fHs. t,) is 

JA 06 J, 1 

the velocity of the middle of segment s^ with respect to the point of 
attachment, O'. 

The point of attachment, O' as discussed above, is placed at the 
beginning of segment s^, at FfsO.t^) * 0. 

The center of mass of the satellite is placed at the end of segment 
3 n ,»t F(s - L,tj). In all of the following, the difference s j +1 - Sj can 
vary with j, allowing for variable integration step size in s, along the 
tether. 

Develop a subroutine to calculate (K(s,t) for any value of s and t. 
Develop a subroutine that can calculate the following quantities, for any 

value s, at time tj, when given K(Sj,t i ), r(Sj tj) ,-|^(Sj ,tj) and 
3 j;j -1 ,2, . . .n, for s between Sj +1 and Sj. 
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For j * 1,2, , n-1 


r(s,t i ) - rCSjt^ - r ^ s j+i) + 


r(3 j+1 ’V " r(s j ,t i ) 

Vi - 5 j 


Ir^-V 


3rt3 j.i, t i ) 

3t 


♦ - 


3r(s J ,t 1 ) 

St 


- S 4 


j>i j 


( S- Sj) 


K(s +1) - K<S,t ) 

K(s,t.) - K(s t,) + =- . (s-s.) 

J J 1 s j*i ‘ s j J 


3K , . K(3, t 1 ) ~ 

3s iS, V “ s-s , 


'j 


a — r(s,t . ) - r(s ,t, ) 

V ^ .., 1 

33 i 3 - S . 


3 2 r(s,t i ) 

3s 2 


SrCs.t^) 3 r (s^tj) 


3 s 


3s 


s-s 


3 3 r(s,t i ) 

3s 2 


3 2 r(s,t i ) S^CSj.tj) 


3s 


3s 


s-s, 


£ ■ c JiVV - 3 <»- 3 j) 


For j • n (s n <sSl): 


(s - s )(56) 


( 57 ) 


( 58 ) 


( 59 ) 


(60) 


( 61 ) 


( 62 ) 


(63) 
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_ r(S n ,t l^ ’ r ^ a n 

r(».t 1 )-r(.„.t l ) *- n ‘ 2 lLl( 3 . s , 

n n-1 


(64) 


3r(3 n ,t i ) 

St ’ at 


at 


^V’V 


S - 3 . 
i n-1 


( 3 -» n ) 


(65) 


_ _ K(3 n ,t.) - K(3 -1,t.) 

ko.V - « Vl , * — - 1 


n n-1 


(s-sj 


( 66 ) 


3 K K(3,t,) - K(s t.) 

|§< 3 ,ti) " « 


(67) 


Sr,. . . r(Sl V ’ r( V V 

3s 3 ’ t i/ gT5 


( 68 ) 


3 2 r 

3 S 2 


3r(s,t fi ) ar(s n t i . 


3 s 


3 s 


s-s 


(69) 


3 3 r (s.tj) 


a 2 r(s,t i ) 3 2 r ( s 1 1 ^ ) 


3 s 




3 - 3 . 


(70) 


3r(s, t . ) 3r(s t.) 

i!! 3 : v ST 5 sH- 


" 5 s 3 t 


3 - 3 . 


(71) 
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mi.: . * 


calculate 


3 2 r(0,tj) 

3t 2 


ff'om eqts (25.1) and (25.2) and calculate 


_ # 

(72 > 


r* C 0 , t ^ ) a) (t^) and u (t^). 

If circular motion is defined, then 


r ■ ak 


(73) 

• 

F * 0 

% 

(74.1) 

• 

ui ■ 0 


(74.2) 

and « is given by Eqt. (33) 




Then, 


with (60), calculate (46) 

- with (56, (57), (64), (58) or (67) calculate (48), 
with (63) or (71) calculate (49), 

with (58) or (66), (57) or (65), calculate (51), 

- with (59) or (67), (60), (57) or (64), (71) or ( 63 ), 
or calculate (52) 

with (60) or (68), (65) or (72) calculate (53) 
with (60) or (68), calculate (54) 
with (62) or (70) calculate (55) 
with (61) calculate (47) 

(end of subroutine) 
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B. Differential Equations 

The differential equations are as follows. From Equation (37) i 


3r(s, t) _ — . 


3t 


v(s,t) 


(75) 


— 2 —*, 

- X(9,t) u(s,t) ♦ AA(s,t) T(s,t) - - 3 r + g(a,t 1 ) (76) 

From equation (35) s 


- v(l,t) 


(77) 


3v(l, t) 

3t 

where, from eqts. (24), (25.1) and (25.2), 


2 -» 


-ll- T(L.t) . + U(L.t) 

m SAT 33 &t 2 


(78) 


3 2 r*(0,t) 3 r s.v. ^ 3 r rel 


3t* 


3t‘ 


at* 


(79) 


m 


J — T(0,t) ♦ —1 


-K 


s.v. 


m s.v. 3 - v - 


(t) 


♦ n 


s.v. 


x (r*(0,t) 



(t)) 


+ a S.V. x ( * 3 .V. X < r * (0 - t ) - r s.V. (t)) 


( 80 ) 
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If Assumption A,B,C,D, and E of Section 5 are made for two-body 
motion, have 


S 2 r*(0.t) 


at* 


_* 

G m Q r (Q,t) 

|F* (0,t)| 3 


And for circular motion with radius a 


3— » 

rr (0,t) 

at 2 


*G m. 


R — with u • 0 and u> given by equation (33). with 


ft* given by equations (31) or (32). Also, with the above 


_» 


assumptions, A,B,C,D and E, the solution for r (0,t) is known analytically 
from the two-body problem. 


If these assumptions are not made, but complete generality is 
retained, eqt. (80) and eqts for Q vec.ff, F must be integrated 

o* V • 

simultaneously with eqts. (75) A (76), to provide 3 2 r — * ^ 

at^ 


From equation (38) 


3T 


* U(5 f t) 


(81) 


3U,„ C(s, t) , H(s, t) ^ . FF( s, t) ^ . F(s, t) , a „v 

^ (3 ’ u ' STvtT * CTTtT TU ' U i(a.t) u(a ' t) * srTTtT (82) 

For the end of the tether, adjacent to the satellite, assuming the 

satellite is a point mass at the end of the satellite, from equation (39), 

letting « T(s»L,t),U s -U(s,t) s 
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C(L,t) + EE( L, t) ♦ -K . 

08 


lira . 


1 T * 1 u .5n 4 t £l 

“sat ‘ • • * p u 3 * a» ^2 


♦ ~~ K*(l,t) - 1 K(s 
“SAT 0 


,tj|. 


Discretizing equation 83 • i for the numerical integratin, and solving 

gives: 


*2- 


- T. P . - ^Lu in . 3f3s - 1 T is . . |n’ 

®SAT ^ p s3s p Sj^Z §3 


--T, 4r • K'(L.t) .iff. ?(».t) 


"s*T L3s 


- (L-s) C (L, t) - (L-s)EE(L.t) 


For the enl of the tether, adjacent to the space vehicle at the 
of attachment, O’, from equation (41) 


C(0,t) + EE(0,t) 



211 1 

•Is^O s 


-h T(s) 


3r(s) 

ds 


*2—7 


- T T(s) — r( ., 8) ♦ 1 K(s,t) 


3s 


3r 3 2 r* (0,t) 

“ 38 ’ at 2 


(83.1) 
for T l 


( 83 . 2 ) 

point 


(84.1) 
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Discretizing equation (84.1) for the numerical integration, and 
solving for U(s) • gives i 




C(0,t) - EE(0, t) * . 1 r±0 j t) m - ( 1 ) . I*(s,t) 


at 


a ' - as * p 


. i( 1) ^lifE^TC) 

p s 3s u *. 


1,-1. 3r(s) 3T(s) 3r(s) 

s 1 p ; * 3s "Is * 3s 


Using eqts. (46) thru (55.1) 


3T(s) . aa^s; . a\s; v/.\ _ r-\ uvu; _ /_\ &«nu; 

~ 5a " ■ — T(,) * mr Ir< ' ,) <s) rsrrr c,) tsw 


AA(s) A(a) 


A(s) 


C(0) 


EE(0) 


{s) as * zm 


3 2 r(0 2 t) 

3t 2 


(84.2) 


In equation (80), define a function 


S(0,ff,K # ) - —2 — K* „ (t) 

m s.v. s ‘ v * 


♦ 2. „ x (r ( 0 , • ) - r m „ (t)) 


s. v. 


S. V. 


* D a.y. * °».v. * (r,(0>t) - r I.v. (t,) 


Then 
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. 5< 0 , S.K-J . _ 

3t 3.V. 


Substituting equation 84.4 into 84.2 gives 


3T(s) 

~5i“ 


. AA (s) . A(s) t>/ _ \ , A( s) t 7 / _ \ 

p gTsT T(s) + ^T K(S) 


,c(0) .EE(O ) 
3 g(st ■ 


♦ s5 ^5tWT s(B -“- k,) 


1 3r(0) 


C Boundary Conditions at s ■ L 


Near the satellite s * 1, we need to impose certain constraints. Let 
the satellite be a point mass, m^,^,, placed at the end of the tether, at 
point s - L. Let s^ denote this point (s^ - L). 

Make the following assumptions. A. Assume that close to the 
satellite, the quantity 


U(s,tj) * (s,t. ) * c(t) 


where c is a constant, with respect to s (but can be a variable in time). 


B. Assume that close to the satellite the quantity ( -v- ) is a 

33 

constant, also with respect to s only. 
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With these constraints, a boundary condition can be derived for the 
tension of the tether at s^. The constraints imposed by A and B, are 
approximate. The constraints are obtained by examination of one 
equilibrium solution near the satellite. Hence when 3*s L but 3-S^ small, 
the error due to constraints A and B should be small when the acceleration 
of the tether and satellite is small and close to an equilibrium solution. 

If the acceleration of becomes significant, than smaller step- 
sizes in s near the satellite should be chosen for the following. 

Let U - tKs.tj), T n - T(s n ,t i ),U n - U(s nj t.),T L - TCs^tj), 

T n-1 * T ^ s n-1 ,t i) U n-1 “ LK s n -i * t i ) and so on, for brevity. The quantities 
T n* U n ,T L’ U,T n-1 ,U n-1 are de P icted in figure 5. 

r 

Applying equation (8/) to segments s n ands n _ 1 gives 


U n _ 1 - - Constant * c 1 (86) 

n n-1 

The quantity T^ is the tension at the point of attachment of the point 
mass m SAT to the tether at T^. 
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For the segment s^ * 


U_ - constant * C, 

" S L ' 3 „ 


Setting (86) and (87) equal gives 


(87) 


T i - T_ T - T , 
L n n n-i 


3, - 3 s - S 
l. n n n-i 


( 88 ) 


And solving for in (88) gives 


3, - S 

X s t + — k Q_ (x _ t 

T n s - s , U n T n-1 ) 
n n-l 


(89) 


Similarly, from Assumption 3: 


X 

SrU.t) 

3s 


3r(s n ,t) 

3s 


3 X ~ 3 n A r(s n » t) ^^n-l* 15 
+ 3 n * 3 n-1 \ 93 


3s 


(90) 


Since the slope, , is assumed constant near s = ^ , other equations 

analogous to equation (89) are : 


T n ■ \- 2 ' 3 < V. ' T n-2> 

n-i n-2 


(91) 


or 


T n- Vi * r Y-~a~; ‘Vi ' T n-2> 

n-1 n-2 


(92) 


The numerical solution procedure is as follows: 


approximated, and 


* 


The numerical integration begins at s - 1 , with T 1 

U 1 given in the following section of this report. The integration of 

equations (81) and (82) procedes to S . Then T. is obtained from equation 

n l 

(83.2), and from (89). 

The value T 1 is adjusted iteratively so that the two values of T n 
coincide to the desired accuracy. 

Equation (90) can be used in two ways. First, the position vector 
along the cable can be initially set so that equation (90) is satisfied. 
Then r(L,t) is allowed to be determined from equations (77) and (78). 

Or Equation (90) can be improve. i rontinually as a constraint on the 
solution for F(L,t) obtained from equations (77), (78). 
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D . Boundary Condition at js-£ 

At s=0, near the point of attachment to the space vehicle, it will be 
assumed that T(s,t) is a constant, or that U(s,t) n - 0. This is the case 
for an equilibrium solution, and we assume that we are close to an 
equilibrium solution. Let thi3 boundary condition be denoted as boundary 
condition II. 


E\ Algorithm for the Solution 

The boundary conditions can be satisfied by various boundary value 
solution methods. We will discuss a standard method often referred to as a 
boundary value solution method using successive solutions of initial value 
problems. Sometimes the method is informally referred to as the "shooting 
method of solution." The first step is to approximate the tension at the 
point of attachment, T(0,t). Then with the boundary condition II, 
U(0,t) * 0, equations (81) and (82) can be integrated from s-0 to s=s 1 
T(0,t) is known and U(0,t)«0 i3 the initial condition for the integration. 

The equations (81), (82) and (83) can be integrated for U(s. t.) and 

J i 

T(Sj, tj), for all s at a fixed time t.. 
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